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Abstract 

Experiments showing a seasonal variation of the nuclear decay rates of a number of different nuclei, and decay anoma- 
lies apparently related to solar flares and solar rotation, have suggested that the Sun may somehow be influencing 
nuclear decay processes. Recently, Cooper searched for such an effect in 238 Pu nuclei contained in the radioisotope 
thermoelectric generators (RTGs) on board the Cassini spacecraft. In this paper we modify and extend Cooper's 
analysis to obtain constraints on anomalous decays of 238 Pu over a wider range of models, but these limits cannot be 
applied to other nuclei if the anomaly is composition-dependent. We also show that it may require very high sensi- 
tivity for terrestrial experiments to discriminate among some models if such a decay anomaly exists, motivating the 
consideration of future spacecraft experiments which would require less precision. 
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1. Introduction 

In a recent paper [1|, early data from a sample of 
137 Cs on board the Messenger spacecraft enroute to 
Mercury were analyzed to set limits on a possible solar 
influence on nuclear decay rates. This work was moti- 
vated by the suggestion put forward in a recent series 
of papers which cite evidence for a drop in the count 
rate of 54 Mn during a solar flare (2), for a correlation 
between decay rates of various isotopes and Earth-Sun 
distance ffl E H [8] [TO] \TT\ \TZ\ , and for period- 
icities in decay-rate data associated with solar rotation 
|fT3l [T4l . Although the suggestion of a solar influence 
on nuclear decay rates has been challenged by the ap- 
parent absence of decay anomalies in some isotopes that 
have been studied lfT5l[l6l[T7ll . and by a recent reactor 
experiment [18|, there is no a priori reason to assume 
that all isotopes should be equally sensitive to a putative 
solar influence, or that the antineutrinos produced in re- 
actors would be the dominant agents through which a 
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solar influence would be exerted. As we have noted 
elsewhere [4, 19, 20 1, the very same properties of de- 
caying nuclides that are responsible for the broad range 
of observed half-lives (e.g., nuclear and atomic wave- 
functions, 2-values, selection rules) would likely render 
nuclides sensitive in different degrees to a putative solar 
influence. 

The suggestion that the Sun is responsible for vari- 
ations in decay rates can be tested directly by studying 
the decay rates of appropriate nuclides located on space- 
craft traveling through the solar system. While such a 
specifically-designed mission has yet to be carried out, a 
number of spacecraft have been launched to date carry- 
ing radioactive nuclides that can be used to constrain de- 
cay anomalies. As mentioned above, Ref. [1J develops a 
general formalism for constraining decay anomalies for 
nuclides placed on board spacecraft, such as the 137 Cs 
source on board Messenger. Additionally, as noted by 
Cooper [16 1, the radioactive nuclides (e.g., 238 Pu) used 
to generate electrical power on spacecraft like Cassini 
via radioisotope thermoelectric generators (RTGs) can 
also be used to set limits on decay anomalies. The goal 
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of this paper is to apply the general formalism devel- 
oped in Ref. [ 1 1 to spacecraft-borne RTGs. 

The outline of our paper is as follows: In Sec. [2] we 
present the general formalism that develops the phe- 
nomenology characterizing a wide range of models of 
solar-induced decay anomalies, and show why experi- 
ments on board spacecraft may be crucial in studying 
them. In Sec. [3] we describe how this phenomenology 
can be applied to obtain limits from spacecraft-borne 
RTGs. In Sec.|4]we apply our formalism to the Cassini 
mission, modifying and extending Cooper's work |[T6l 
to a wider class of models. We conclude in Sec. [5] with 
a summary and discussion of our results. 



2. General Formalism 

We begin by briefly reviewing the formalism devel- 
oped in Refs. |[T] |4) to describe anomalous radioactive 
decays. If N(t) is the number of unstable nuclei in a 
sample located at position r, then we will assume that 
the total activity of the sample is 



dN(t) 
dt 



= A(t)N(t) = [A + A 1 (?,t)]N(t), (1) 



where A$ represents the intrinsic contribution to the de- 
cay rate of the unstable nuclei, along with a possible 
time-independent background arising from new interac- 
tions. Here A\{r, t) <K Aq characterizes the anomalous 
position- and time-dependent contribution to the decay 
rate, assumed, in this case, to arise from the Sun. 

As noted above, there is evidence using terrestrial ra- 
dioactive samples suggesting that decay rates are corre- 
lated with sample-Sun separation. To study this effect, 
we will assume a specific phenomenological form for 
Ai{?,t) given by |Q] 



A x (?,i) = Ao? 



■(») 



R 



r(t) 



(2) 



where n = 1,2, 3, . . ., r = |/*|is now the sample-Sun sep- 
aration distance, R = 1 AU = 1.495979 x 10 8 km, and 
£^ n) is a composition-dependent dimensionless parame- 
ter characterizing the "strength" of the decay anomaly 
for a specific nucleus. The form of Eq. (|2]i is designed 
to encompass a broad range of theories. If the decay 
anomalies are caused by a flux of neutrinos from the 
Sun or by an inverse-square law field, then one expects 
n — 2. On the other hand, radioactive decay processes 
might be affected by one of the many proposed new 
long-range inverse-power-law interactions. In this case, 
the potential energy between point particles of mass m\ 



and ni2 separated by distance r may be written in the 
general form ETI 



(3) 



where a„ is a dimensionless constant, G is the New- 
tonian gravitation constant, and ro = 10~ 15 m is a 
length scale chosen by convention. For example, n — 2 
can arise from the simultaneous exchange of 2 mass- 
less scalar particles 11221 . n — 3 from 2-massless- 
pseudoscalar exchange ||231 l24l . and n — 5 from 2- 
neutrino [25, 26 1 and 2-axion exchanges [24|. 

Because the Earth-Sun variation of r is small, it is 
difficult to obtain the dependence on the power n from 
a purely terrestrial experiment, which is why the use of 
spacecraft with widely varying values of r is important. 
To show this, we substitute Eq. Q into Eq. ([T]), giving 



dN(t) 
dt 



R 



N(t). (4) 



This form is not very useful since it depends on the in- 
stantaneous value of N(t). Following Ref. (H, we can 
eliminate this dependence by first integrating Eq. Q, 
obtaining 



N(t) = N exp {-A [t + £ n) I (n \t)]} , 
where No = N(Q) and 



f dt' 


R 


Jo 


r(t') 



I m (t) = 

Then differentiating Eq. |5]) gives 
- N(t) = AoN expMo?) 1 1 + f 00 



(5) 



(6) 



R" 



A I {n) (t) 



(7) 

where N(t) = dN(t)/dt and assuming g n) At <K 1,. For 
our application, we will set t — to be the launch time 
of the spacecraft, in which case r(0) = R. Then, using 
J ( "'(0) = 0, Eq. Q yields 



-#(0) = AoiVb[l +£<*>]. 



(8) 



For practical purposes, it is useful to express Eq. |7]l 
in terms of directly measurable quantities. If we define 
A = Aq[\ + to be the decay rate observed on Earth 
when r = R, then ^ = A[l - % (n) ] + 0[^>f. Substitut- 
ing this into Eq. (j7]) and retaining only terms of leading 
order in £ n \ we find after using Eq. <JsJ that Q]] 



N(t) N(0)e- /h [l + % {n) B { "\t)\ , 



(9) 



2 



where 



3. Application to RTGs 



B {n \t) = 



R 

W) 



l\-A[l ( "\t)-t\. (10) 



For later purposes, we note from Eq. ( 10 1 that B ( "\Q) 



0. The terms in braces in Eq. ( 10 1 represent the anoma 



lous contributions to N(t) arising from the circumstance 
that at some t the sample is at r(f), rather than at R. The 
second term in square brackets represents an additional 
cumulative contribution to N(t) from the variation of 
r(f) from t = to t, and this is generally nonzero even 
when r(t) = R (see below). 

We will now demonstrate why spacecraft missions 
may be needed to distinguish among the inverse-power- 
law models given by Eq. Q. Since the eccentricity of 
the Earth's orbit is small (s B 0.0167), the Earth-Sun 
separation can be written as 



■Jt) = + sin(^| + {Ws. (II) 



where r ffl R is the mean Earth-Sun separation, T = 1 
year is the orbital period, and t = when the Earth is at 
r — R. Substituting Eq. (Hi into Eq. (|6]( and integrating, 
while keeping only terms of 0(s @ ), we find 



ns e T 
2?r 



cos 



2nt 



(12) 



so B^"\t) given by Eq. ( 10 1 becomes, after neglecting 
terms of O(sl), 



B ®\i) - ~ ns <s \ sin 



2nt\ AT 
~2~k 



1 - cos 



2nt 



(13) 

Combining Eqs. (|9| and ( 13 ), we find that the activity 
of a sample on Earth is given by 



N m (t) - N @ (0)e- M \l-n^e @ isin 



2nt 



AT 



1 - cos 



2nt 



(14) 



Since n and £^ appear together only in the combination 
nt*- n) in Eq. ( 14 ), it is not possible to distinguish between 



theories with different values of n using terrestrial decay 
experiments unless one has sufficient sensitivity to de- 
tect effects of 0{s 2 @ ) ~ 3 x 10 4 . On the other hand, by 
placing radioactive samples on board spacecraft, a de- 
cay anomaly can be probed over a wide range of r(t), 
allowing one to discriminate more easily between the 
various powers of n given in Eq. |2]). 



The thermoelectric generators on board spacecraft 
use the Seebeck effect to convert heat from a hot reser- 
voir into electrical power, while waste heat is exhausted 
into a cold reservoir. In a simple model lF27l . the total 
electrical power output can be written as 



P d (0 = ZAV tt 



(15) 



where / is the current, and the terminal voltage across 
the generator is given by 



AVterm = ce n , p (T H - T c ) - IRh 



(16) 



(17) 



Here a, up is the differential Seebeck coefficient for the 
n- and /?-doped semiconductor generator legs, Th and 
Tc are the temperatures of the hot and cold heat reser- 
voirs, and R m is the internal resistance of the generator. 
The efficiency of a thermoelectric generator is usually 
expressed as 

_ p * 
riTG - -r-, 
Qh 

where Qh is the thermal power input from the hot reser- 
voir. 

An RTG is a thermoelectric generator that uses a 
radioactive material as its hot temperature reservoir. 
RTGs have been placed on board spacecraft since the 
early 1960s lF27l l28l . and are typically used on probes 
traveling outward in the solar system where solar panels 
do not provide sufficient power. These missions include 
Pioneer 10 and 11, Voyager 1 and 2, Galileo, Ulysses, 
Cassini, and New Horizons. In addition, RTGs were 
used to power instruments on the two Viking landers on 
Mars, and instruments on the Moon placed by Apollo 
astronauts. 

The thermal power P±(t) generated by the radioactive 
heat source of an RTG is directly proportional to N(t), 
the activity of the radioactive material used. Therefore, 
we can use Eq. ([9]) to relate the thermal power generated 
at time t to the thermal power produced at launch t = 0: 

AhW = P t h(0)e- At [l + t (n) B {n) (t)\ . (18) 

Unfortunately for our purposes, this thermal power is 
not directly observed. Instead, the electrical power 
f e i(0 given by Eq. ( [T5] ) is measured. To relate P e i(0 
to Pth(t), we follow Cooper ifTBI and introduce the di- 
mensionless efficiency function e(f) defined by 



where 



e(t) = T]TG(t)lmd(t), 



(19) 



(20) 
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and 



Qh 



(21) 



is the fraction of the total radioactive thermal power that 
flows into the generator. Eq. ( 18 1 can then be rewritten 

as 

Pd.it) = P d (0) 6(0 e~* [l + t {n) B {n) (t)\ , (22a) 
= P e i(0) e(t) 2-' /r " 2 [l + Z in) B in) (t)\ , (22b) 

where t\/2 = (ln2)//l is the half-life of the nuclide. We 
note that e{t) is time dependent since RTG efficiency 
generally decreases with time. 

For our problem, eit) must be determined in a man- 
ner that is not influenced by the presence of a decay 
anomaly, and one option is to use a computer model. 
However, the model commonly used to describe RTG 
performance (DEGRA 11301 ) assumes the usual expo- 
nential decay of the radionuclide, preventing its use in 
testing the exponential decay hypothesis (16). The al- 
ternative suggested by Cooper ifTBI is to use an empiri- 
cal approach that capitalizes on the fact that there may 
be points along a spacecraft's trajectory where anoma- 
lous decay effects make no contribution, which from 
Eq. (22 1, occurs whenever B ( "\t) = 0. In his analy- 
sis of Cassini power data, Cooper utilized five points at 
times where r(t) = R, but from Eq. (JTOj we see that 



B [n \r = R)=-A [l (n \t) - t\ , 



(23) 



which does not generally vanish for t > 0, though in 
the specific case of Cassini studied by Cooper, this is 
a reasonable approximation. If we define the decay- 
normalized electrical power as 



Pit) = 2 +,/ '" 2 



PdO) 



Pdt) 

Pel(0)' 



(24) 



and use only power data where B (/,) (f) = 0, then we can 



solve Eq. (22 1 for e{t) in terms of known quantities, 

eit) = Pit; B (n) = 0), (25) 

and fit the results to an empirical function as we will 
demonstrate in the next section. 

This analysis indicates that decay anomalies can, in 
principle, be detected using the electrical power output 
as a proxy for the (time-dependent) energy release from 
nuclear decays in spacecraft RTGs. All that is required 
is: (1) a spacecraft trajectory with a significant varia- 
tion of r(t) (preferably including r{t) <K R), (2) accu- 
rate measurements of the spacecraft's position and elec- 
trical power production, and (3) a good model of the 



RTG efficiency function e(t) that does not assume the 
radioactive decay law. If a computer model for the ef- 
ficiency is unavailable, an empirical approach can be 
used, provided there are a sufficient number of points 
where Bit) = 0. The original Solar Probe Plus mis- 
sion, which would have used a spacecraft powered by 
RTGs while using a Jupiter gravity assist maneuver to 
place it into solar polar orbit, would have been ideal for 
our purposes BP . However, this planned mission has 
been changed to avoid the use of RTGs. Of the inter- 
planetary missions that have used RTGs, only Cassini 
and Galileo actually crossed the 1-AU orbit radius af- 
ter launch, allowing the empirical approach. Since the 
Cassini mission crossed the most times, it provides the 
most complete data to model the RTG efficiency empir- 
ically and so will be used in this paper. 

We note in passing that accurate modeling of RTGs 
has also been important in understanding the Pioneer 
anomaly, the small anomalous acceleration of both Pi- 
oneer 10 and 11 spacecraft [32] • m order to explain 
the observed temporal decay of this acceleration, the 
changes in thermal recoil forces on these spacecraft due 
to the degradation of the RTGs need to be accurately 
modeled IT331. 



4. Application to Cassini RTG Data 

Launched in 1997, the Cassini spacecraft made sev- 
eral gravity-assisted flybys of Venus, Earth, and Jupiter 
before entering the Saturn system in 2004. In principle, 
this nearly 7 -year-journey provides both a sufficiently 
long duration and a substantial variation in r(f) to ren- 
der Cassini ideal for a test of solar-influenced radioac- 
tive decays. However, following Cooper IfTBI . we will 
only use the first two years of data, when 0.6732 AU 
< r < 1.6215 AU, because only within this range can 
we determine with some confidence the efficiency func- 
tion eit). For t > 2 years, there are no longer any times 
where B ( "'(f) = which can be used to perform the em- 
pirical fit to determine eit). 

The Cassini spacecraft uses three General Purpose 
Heat Source Radioisotope Thermal Generators (GPHS- 
RTGs), each producing nearly 300 W of electrical 
power from 572 thermoelectric "unicouples." A GPHS 
generates approximately 4410 W of thermal power at 
the beginning of life from PUO2 fuel pellets enriched 
to about 80% 238 Pu, corresponding to about 8.1 kg of 
238 Pu per generator |34|. The pellets were formed from 
PUO2 powder obtained from Russia, where the pluto- 
nium was created by irradiating 237 Np in a high-flux re- 
actor to form 238 Np which has a half-life of 2.4 days 
(29). A 238 Np nucleus decays to 238 Pu via a B~ -decay. 
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Figure 1 : Electrical power output from each of the three Cassini RTGs 
for the first two years of the mission. The digitization arises from the 
current measurements. 



Figure 2: B*"'(r) of Cassini obtained using Eq. jlO| for the first two 
years of the mission. The five points where B l "\t) = 0, which occur 
when r = R, were used to obtain the empirical efficiency function as 
shown in Fig. [5] 



For the rest of this analysis, we will assume that all 
of the thermal heat production of a GPHS results from 
238 Pu, which has a half-life f;/2 = 87.7 years. (The other 
most likely nuclide present would be 239 Pu, whose half- 
life is 2.41 x 10 10 yr, and hence decays too slowly to 
contribute significantly to the energy production in an 
RTG.) 

Our analysis begins with Cassini's trajectory and 
RTG current 7(f) data provided by the Jet Propulsion 
Laboratory. Since the current is monitored separately 
for each RTG, each can provide a separate determina- 



tion of the material-dependent parameter ^ for Z38 Pu, 
and together the RTGs can be used to assess the valid- 
ity of the empirical efficiency function e emp (f)- A to- 
tal of 157,465 current measurements (in increments of 
0.03959 A) were taken over the 2-year period at irreg- 
ular intervals, and a small number of obviously spu- 
rious points were removed before our analysis began. 
The Cassini power system is regulated with a variable 
shunt radiator to maintain a constant terminal voltage 
AVterm = 30.0 ± 0.2 V across each RTG to maximize 
power production ||35l , and hence the RTG electrical 
power was obtained from the current measurements us- 
ing Eq. (13 1. Since the variable shunt regulator is lo- 



cated "downstream" from the RTG, its operation does 
not affect our ability to determine the RTG's power out- 
put as a function of time. The results are shown in Fig.[T] 
with digitization of the current measurements clearly 
evident. The power decreases more rapidly than ex- 
pected from the radioactivity exponential decay law for 
238 Pu due to unicouple degradation among other factors. 

Using the Cassini trajectory data, B (n) (f) for n = 1-5 
was computed using Eq. ( 10 1 for < t < 2 years, and is 



plotted in Fig. [2] We note that B <n \t) increases with n, 
so the constraints on R£ become more stringent as n in- 
creases. From Eq. (JTOj, we see that the largest values of 
B (,!) (f), which lead to the largest decay anomalies, occur 
when r{t) is smallest. For Cassini, this occurs at its clos- 
est approach to the Sun, r = 0.6732 AU at t =s 0.45 yr, 
with smaller peaks at f ^ 1.7 yr. We also note from 
Fig.gthat there are five points where B ( "\t) = 0, which 
will be needed to obtain the empirical efficiency func- 
tion. 

Using the data from each of the three RTGs for 
the first two years, we calculated the decay-normalized 



power !P(f) using Eq. (24i. For the 5 data points where 
B ( - n \t) = 0, we then determined e(f) following Eq. (25 1 
by plotting for each RTG e(f) = P(t; B = 0) versus f, 
and then fitting the results to an empirical function sim- 
ilar to that suggested by Cooper ifTBI . 



c emp 



(f) = 2~' /7; ' lr (l + at + bt 2 + ct 3 ) . 



(26) 



where T e g, a, b, and c are constants. [Note that our 
procedure differs from Cooper who only used e emp (f) = 
2~'' T " S for the points where r(f) = R. The third-order 



polynomial portion of Eq. ( 26 1 was introduced later and 
fit using all of the trajectory data < t < 2 yr.] The 
result for n — 2 is shown in Fig. |3j the results for 
the other values of n are virtually identical. We see 
that there is some relative variation among the three 
RTGs, and we use this variation to estimate the rela- 
tive uncertainty in the empirical efficiency function to 

be (56emp/^emp 

< 0.5%. 

We can assess the validity of our empirical approach 
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t [yr] 



Figure 3: Fits of Pit) where = to e e mp(0 given by Eq. |26) 

for the three Cassini RTGs for n = 2. A small variation in efficiency 
among the RTGs is observed. 



for determining e(t) by examining the quantity 

~ nt) 



A(f) = 



£emp(0 



- 1 



(27) 



which would vanish if e e mp(0 gave a perfect character- 
ization of P(t). In Fig. H] A is plotted versus t and r. 
Using the combined data from all three RTGs, we find 



|A(0I < 0.005, 



(28) 



with no significant systematic dependence of A on t or r 
other than the effects due to the digitization of the power 
seen in Fig. [T] Thus, e emp (f) gives a good characteriza- 
tion of the decay-normalized power of the RTGs to the 
level of 0.5%. Systematic effects due to solar heating 
of the RTGs, which would depend on 1/r 2 , are not ev- 
ident. The semiconductor legs of the unicouples gener- 
ating the power via the temperature difference Th - Tc 
are only 20.3 mm long 11341 . Only a differential heating 
over this distance scale would produce a change in the 
RTG electrical power production, which is not seen. 

Now that the efficiency function e(f) = e emp (f) has 
been determined, all quantities in Eq. (221 are known 



except for Inserting Eq. (24 1 into Eq. ( 22 1 and set- 
ting e(f) ^ e emp (t), we can solve for to obtain 



An) 



A(rt 



(29) 



where we have used Eq. ( |27]i . 
limit on A(f) given by Eq. (|28|> 



Inserting the numerical 



into Eq. ( 29 1, we find 



< 



0.005 



(30) 



0.005,) 



W!?"ffff, 




t(yr) 
(a) 



RTG 1 
RTG 2 
RTG 3 




1.2 
r(AU) 



(b) 



Figure 4: A given by Eq. \21) plotted versus t (a) and r (b). The 
striped patterns are due to the digitization of the power seen in Fig.^ 
No other significant systematic pattern is observed. 
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Table 1 : The constraints on for Pu obtained from Eq. |30| when 



n 


R («) 

D max 






1 


0.48 


< 1.0 x 10" 


-2 


2 


1.19 


< 4.2 x 10" 


-3 


3 


2.25 


< 2.2 x 10- 


-3 


4 


3.82 


< 1.3 x 10- 


-3 


5 


6.14 


< 8.1 x 10- 


-4 



The most stringent limits on £r , which are obtained 

when |-S (n) (f)l = #max» me maximum observed value of 
B {n \t), are presented in Table [I] We see that the tight- 
est constraints are obtained for the largest values of n, 
which is consistent with our earlier discussion, while a 
comparitively poor constraint is obtained for n = 1 since 
i?max is relatively small. 



5. Discussion 

From our analysis of the Cassini RTG power data, we 
find no evidence of a solar-influenced decay anomaly 
for 238 Pu having the inverse-power-law form given by 
Eq. (j2j) for n = 1-5. Unfortunately, it is difficult to com- 
pare our constraints on ^ given in Table [T] with the 
earlier analysis by by Cooper [16|, who also found no 
anomaly. Schematically, Cooper fit the total normalized 
RTG power production to an empirical efficiency func- 
tion, and then fit the resulting residuals to two different 
possible functions characterizing the decay anomaly: 



residuals = 



P 



(31) 



where R = 1 AU. This procedure yielded the 90% con- 
fidence level limits 



\a\ < 0.84 x 10" 
LSI < 0.99 x 10" 



(32a) 
(32b) 



However, despite their superficial similarity, Cooper's 
formulas given by Eq. 



(31 



are not directly related to 
our more fundamental formula for Ai given by Eq. (|2]). 
The anomaly in the power production of the RTGs in 
our approach is not a simple power law since one needs 
to take into account the fact that the sample's activity 
at time t actually depends on the sample's entire history 



because of its dependence on the function P"\t) given 
by Eq. ([SJ. Since the two approaches of characterizing 
the power anomaly are so different, it is difficult to relate 



Cooper's limits on a and B given by Eq. ( 32 1 to our con- 
straints on ^ and respectively, given in Tablejl] 

As noted in the Introduction, one expects an anoma- 
lous decay mechanism to depend on the nuclei and de- 
cay process so one cannot, without additional assump- 
tions, use the limits on £, {n) for 238 Pu to constrain anoma- 
lous decays of other nuclides or to refute the obser- 
vations of previous experiments. The Cassini RTGs 
are powered exclusively by alpha-decays due to the ex- 
tremely long half-lives of the uranium daughters, partic- 
ularly 234 U (f 1/2 = 246, 000 yr) which is the first daugh- 
ter of 238 Pu. The absence of a decay anomaly for 238 Pu 
contrasts with the other isotopes in which an anomaly 
has potentially been observed Gl |4] [3 El O [HI EJ UHl HH 
fl2l [131 f!4l . that are beta-decays or (like 226 Ra) where 
one actually measures a significant beta-decay compo- 
nent of daughter products |[36l . For example, Parkho- 
mov |S), found time-dependent fluctuations in 90 Sr, 90 Y, 
and 60 Co, all of which are beta decays, but not in 239 Pu 
which, like 238 Pu studied in this paper, is an alpha de- 
cay. Each radioactive nuclide needs to be examined sep- 
arately for anomalous decay processes, and £™ deter- 
mined for each. If Eq. (|2]) is correct, all anomalous ef- 
fects should be characterized by the same power n, and 
all experiments using the same nuclide should yield the 
same value of Work is currently under way to ap- 
ply Eq. Q to the results of terrestrial experiments, as 
described in Sec. 2. The only previously reported result 
is = (2.8 ± 8.1) x 10~ 3 for the 137 Cs sample aboard 
the Messenger spacecraft [T]. 

We have also shown that if an inverse-power-law 
form given by Eq. Q exists, terrestrial decay experi- 
ments will require unusually high sensitivity to discrim- 
inate £ (n) from n since only the combination nt^ n) ap- 
pears to first order in the orbital eccentricity. Thus, 
experiments conducted on board spacecraft may be 
needed to distinguish among the various powers of n. 
While we and Cooper ifTBI have shown that spacecraft- 
borne RTGs can be used, the isotopes used for power 
generation are very limited, and hence dedicated exper- 
iments using nuclides that have demonstrated a decay 
anomaly in terrestrial experiments should be used. 

Since RTGs will continue to be a power source on 
board spacecraft, improved instrumentation on future 
missions could not only be used to refine RTG perfor- 
mance models, but also provide improved tests of the 
radioactivity exponential decay law. The two major fac- 
tors limiting constraints on decay anomalies using the 
Cassini RTGs are the uncertainties of the empirical ef- 
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ficiency model and the resolution of the electric current 
(and hence, power) measurements. Additional thermal 
measurements would also be useful in allowing a more 
direct determination of the thermal power output of the 
RTG heat source. 
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